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Abstract. The B-meson decay constant /b has been calculated from unquenched lattice QCD in the 
unphysical region. For extrapolating the lattice data to the physical region, we propose a phenomenological 
functional form based on the effective chiral perturbation theory for heavy mesons, which respects both the 
heavy quark symmetry and the chiral symmetry, and the non-relativistic constituent quark model which is 
valid at large pion masses. The inclusion of pion loop corrections leads to nonanalytic contributions to /s 
when the pion mass is small. The finite-range regularization technique is employed for the resummation 
of higher order terms of the chiral expansion. We also take into account the finite volume effects in lattice 
simulations. The dependence on the parameters and other uncertainties in our model are discussed. 
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1 Introduction 

An accurate determination of the Cabibbo - Kobayashi - 
Maskawa (CKM) matrix of the Standard Model and tests 
of its consistency and unitarity constitute an important 
part of current research in particle physics. Among the 
most important matrix elements is the _B-meson decay 
constant /s, which is needed to determine the CKM ma- 
trix elements such as Vtd- Lattice Quantum Chromody- 
namics (QCD) simulations, which provide a way to de- 
termine fs from the first principles of QCD, are of fun- 
damental importance since the determination of the B- 
meson decay constant still remains beyond the reach of 
experiment. 

In Ref.[lJ, the authors studied fs extensively from full 
and partially quenched lattice QCD. In the simulations 
they applied two actions, the unquenched gauge configu- 
rations by the MILC Collaboration[2 and the improved 
staggered light quark action [3], for simulating the light 
quarks {u,d,s). The non-relativistic QCD (NRQCD) for- 
malism was used to treat the heavy quarks (c, b) [U . The 
lattice data were obtained in the region where the mass 
of Ught quark is about in the region 20 MeV ~ 80 MeV[I] 
(corresponding to the pion mass between around 250 MeV 
and 500 MeV), which is larger than its physical value. The 
light quark masses in these simulations are much closer 
to the physical region than previous work. This makes it 
possible to reduce the uncertainty from chiral extrapola- 
tion. To extrapolate the lattice data to the physical region 
the staggered chiral perturbation theory was employed in 
Refs.dld. 



In QCD, when the heavy quark mass goes to infinity, 
the strong interaction can be described by the heavy quark 
effective theory (HQET) [6] , which contains the heavy quark 
symmetry SU{2)f x SU{2)s. On the contrary, when the 
masses of light quarks (w, d, s) approach zero, the strong 
interaction has SU{3)l x SU{3)u chiral symmetry which 
is spontaneously broken to the SU{3)v subgroup, leading 
to eight pseudoscalar Goldstone bosons. When the light 
quark mass (or the pion mass) is small enough, the effec- 
tive chiral Lagrangian[7 containing the chiral and heavy 
quark symmetries may be applied, based on power count- 
ing which is the foundation of chiral perturbation theory. 
The authors of Ref. W found that the regime where power 
counting may be applied is < ttItt < 180 MeV by analyz- 
ing the chiral extrapolation of the nucleon mass to fourth- 
order in the expansion and demanding the accuracy of the 
prediction for the nucleon mass to be one percent. This 
regime is beyond the reach of current lattice simulations. 
In order to extrapolate lattice data which are outside the 
power-counting regime, the so-called finite-range regular- 
ization (FRR) technique is proposed [SlI MTUUTT] . in which 
a finite-range cutoff (which is physically of order of the in- 
verse size of the pion source) in the momentum integrals of 
the pion loop diagrams is introduced for the resummation 
of the chiral expansion. It has been shown strictly that 
FRR is mathematically equivalent to minimal subtraction 
schemes such as dimensional regularization to any finite 
order[8]. Furthermore, the model dependence associated 
with the shape of the regulator is at one percent level in 
the range < < 1 GeV, which is well outside the 
power-counting regime. In the present work we will adopt 
FRR while applying the effective chiral Lagrangian to ex- 
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trapolate the lattice data for the _B-meson decay constant 
since these data are outside the power-counting regime. 

The constituent quark model has been shown to work 
quite well although it is a very simple model where Quan- 
tum Chromodynamics is not employed. At large pion masses, 
we employ the non-relativistic constituent quark model to 
fit the lattice data. In fact, for all the hadron properties 
(such as heavy meson masses and decay constants) which 
have been calculated in lattice QCD, the lattice data vary 
slowly and smoothly when the light quark mass is larger 
than 60 MeV or so. This is characteristic of constituent 
quark behavior and suggests that in this region the con- 
stituent quark model should be most appropriate [l2]. In 
other words, the hadron properties vary as a function of 
the constituent quark mass in the region of heavier quark 
mass. It has been shown that the constituent quark model 
is consistent with the modeling of QCD via the Dyson- 
Schwinger Equation in ladder-rainbow truncation when 
the quark mass is large[TT]. Obviously the lattice data 
at large light quark masses in Ref.|lJ appear in the re- 
gion where the constituent quark model is suitable to be 
employed. Consequently, we apply the constituent quark 
model in this large quark mass region and then analyti- 
cally continue the hadron properties to the physical mass 
regime in a manner consistent with chiral symmetry, be- 
cause of which nonanalytic terms are involved. 

In Ref.[Tj, the authors extrapolated the lattice data 
for the quantity <Pb {'^b = Ib^/ttib) with the following 
formula: 



'b = co{1 + a + analytic terms), (1) 



where Cq is a constant associated with the axial- vector cur- 
rent which destroys B meson, A is the term encompassing 
the chiral logarithms, and the analytic term is linear in the 
light quark mass. It is obvious that when Eq.(IT]) is used 
for extrapolating the lattice data, power counting in chiral 
perturbation theory has been applied when the light quark 
mass (or the pion mass) is small and the term linear in the 
light quark mass has been assumed when the light quark 
mass is large. However, as mentioned before, the lattice 
data appear in the light quark mass range 20 MeV ~ 80 
MeV (or the pion mass range 250 MeV ~ 500 MeV). This 
is already outside the regime where power counting can be 
applied (0 < to^ < 180) for small pion masses and inside 
the regime where the constituent quark model is most ap- 
propriate for large pion masses. Therefore, based on our 
above argument, instead of using the form like Eq.([T]) to 
extrapolate the lattice data for the B-meson decay con- 
stant, we adopt FRR when applying chiral Lagrangian 
and the constituent quark model at small and large pion 
masses, respectively. 

The outline of this paper is as follows. We will give a 
brief review for the chiral perturbation theory for heavy 
mesons in Sect. 2. In Sect. 3, we will propose a functional 
form for extrapolating the lattice data to the physical re- 
gion based on the effective chiral perturbation theory for 
heavy mesons and the constraints from the constituent 
quark model at large pion masses. Then in Sect. 4, we will 
present our fits to the lattice data and give the numerical 



results. Finally, a summary and discussion will be given 
in Sect. 5. 



2 Chiral perturbation theory for heavy mesons 

The chiral perturbation theory for heavy mesons, which 
describes the interactions of mesons containing a single 
heavy quark with light pseudoscalar bosons, contains both 
the heavy quark symmetry SU{2)f x SU{2)s and the chi- 
ral symmetry SU{3)l x SU{3)r, provided mg — > oo and 
ruq ^ {Q = c or b, q ~ u,d, s ), respectively [3. 

The pseudoscalar Goldstone bosons are incorporated 
in a 3 X 3 unitary matrix 



E = exp 



/2iM 
\17 



(2) 



where fj^ is the pion decay constant, ~ 131 MeV, and 



M 



V 



%/2 



7r+ 




K+ 















(3) 



While discussing the strong interaction between heavy 
mesons and light pseudoscalar mesons it is convenient to 
introduce 

f = v^. (4) 
Under a chiral SU{3)l x SU{3)u transformation. 



L^C/t = U£,R\ 

where the unitary matrix [/ is a complicated nonlinear 
function of L, R and the pseudoscalar Goldstone boson 
fields. 

In order to describe the interactions of the Goldstone 
bosons with the heavy mesons containing Qq"^ (here a — 
1,2,3 for u,d,s quarks, respectively), it is convenient to 
introduce a 4 x 4 matrix Ha given in Ref.[7], 



1 + / 



(5) 



where Pa and P*^ are the field operators that destroy a 
heavy pseudoscalar meson {Pa) and a heavy vector meson 
(P*) with four- velocity w, respectively. P*^ satisfies the 
constraint as the following: 



(6) 



Under SU{'i)L x SU{?>)r transformation. 

Ha - H,UI, (7) 
and under the heavy quark spin transformation. 

Ha ^ SHa. (8) 
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Defining 
we have 



15) 



i + i' 



It is convenient to introduce a vector field V^^^ 



and an axial- vector field A^^^, 



(9) 
(10) 

(11) 

(12) 



to establish the Lagrangian for the chiral perturbation the- 
ory. 

The most general form of the effective Lagrangian den- 
sity, which is invariant under the heavy quark symmetry 
and the chiral symmetry and should be invariant under 
Lorentz and parity transformations as well, is as follows^: 

C = -TY[Haiv^{D^H)a] + .gTr(5„i/b7^A^^75), (13) 

where g is the coupling constant describing the interac- 
tion between heavy mesons and Goldstone bosons, which 
contains information about the interaction at the quark 
level. So it cannot be fixed from the chiral perturbation 
theory for heavy mesons, but should be determined by ex- 
periments. The covariant derivative in Eq.(13) is defined 
as 



ba ' 



(14) 



Then we have the following explicit form for the in- 
teraction of heavy mesons with Goldstone bosons after 
substituting Eqs.(l5l)(Il0l)(Ill|)(IT2l) into Eq.(Il3l): 

= j^v^ [M, d^M]UP*alPr - PlPb) 

-^{p:lPbd^M,a + pip:,,d'^M,a 

+,,M-p.p*tp^;„^a^M,,), (15) 

where 0{M^) terms are ignored. 

Taking the mass difference between Pa and P* into 
account, the following term has to be added into C in 
Eq.(13): 



— Tri/.d^'-i/a^Tu., 



(16) 



where A2 is a constant containing interaction information 
at the quark level. 

So the propagators for heavy pseudoscalar and vector 
mesons are 

(17) 



and 



(18) 



respectively, where p is the residual momentum of the 
heavy meson. In Eqs.(17)(18) 



A 



8A2 



(19) 



which is the mass difference between vector and pseu- 
doscalar heavy mesons. 



3 Formulas for the extrapolation of the 
5-meson decay constant 

The decay constant of a pseudoscalar heavy meson, P, is 
defined by 

(0| J^|P) = tfpmpv'', (20) 
where the axial-vector current is 



J'' = g7^(l-75)g, 



(21) 



q denotes u, d light quarks and Q represents heavy quarks 
c, b. 

The current in Eq.(21) can be written in the low energy 
chiral theory as[l3j 



^^' = (f )Tr[7'^(l-75)ff4j, 



(22) 



where a is a parameter. 

Substituting Eq.(5) into Eq.(22), we have the following 
explicit form for the current J^: 



J^ = ^a{P;''-v^^Pt)^l. 



(23) 



We make a Taylor expansion for M and omit 0{-p) 
terms. This leads to the following expression for J^: 

1 



(24) 



The diagrams for the heavy meson decay constant to 
one pion loop order from Eqs.(15)(24) are shown in Fig. 
1. We can see that there are three diagrams for one pion 
loop corrections to the pseudoscalar heavy meson decay 
constant. 

From Fig. 1 (a) , it is easy to obtain 



fi") I 



(25) 



where fp^ is the heavy pseudoscalar meson decay constant 
at the tree level where pion loop corrections are not taken 
into account. 

Fig. 1(b) represents the pion loop correction from the 
axial current itself. It can be expressed as 



Ab) 

J p 



d^fc 



1 



(27r)4 /c2 _ ^2 ' 



(26) 



where k is the four momentum of the pion in the loop 
and is the pion mass which is not necessary to be its 
physical mass. 
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Choosing the contour in the lower half complex plane 
of fco, in which there is only one pole for ka, we have the 
following result after integrating over fco^ 



(0) 



Ab) -^r 



d^fc 1 

(27r)3 wk ' 



(27) 



where Wk = -v/lfeT^ + 

Fig. 1(c) vanishes because of the identity 

v^[9tiu - v^v^) = 0. 

The contribution of Fig. 1 (d) to the matrix element in 
Eq.® is 

9»g^ .(0) y 

multiplied by impv^, where the self-energy contribution 
is proportional to 



3 



d^fc 



(27r)4 [y.[p-k)- iA](fc 



(29) 



Since u -p is a Lorentz scalar, we are free to choose the 
special frame in which the heavy meson is at rest: u° = 
l,v = 0. Choosing the same contour in the evaluation 
of the integral as before, we have (details are given in 
Ref.[10]) 

d^fc IkP 



X = 



6 



(30) 



(27r)3 Wk{wk - 5) ' 

where 5 — v ■ p — jA. 

The contribution to fp from the one-loop wave func- 
tion renormalization, /p'^^ can be obtained from the deriva- 
tive of X with respect to v ■ p: 

1 dX . 



2dvp"'P^ 



Then we have 



Ad) 
Jp 



35 



"■^ AO) 

J (2^)3 u,fe («;,,- 5)2- 



(31) 



(32) 



Adding all the contributions from the diagrams in Fig. 
1, we have the following expression for fp: 



Jp — Jp + Jp + Jp , 



(33) 



where the first term represents the contribution at the 
tree level, and the other two come from one pion loop 
corrections. 

When the light quark mass is near the chiral limit, pion 
loop corrections give the dominate contributions to fp; on 
the contrary, pion loop contributions vanish in the limit 
m,r ^ oo. In order to extrapolate the lattice data, we also 
have to know the behavior of fp at the large pion mass. 
Based on the non-relativistic constituent quark model, it is 
known for a long time that, up to logarithmic corrections, 
fp obeys the asymptotic scaling law when the constitute 
quark masses are very large: 



constant. 



(34) 



The above behavior is also obtained in the Dyson- Schwinger 
Equation approach 14j. 

Therefore, it is convenient to extrapolate fpy/mp in- 
stead of fp. Then we only need to introduce a constant 
parameter while fitting the lattice data in the region where 

is large. 

Based on the above argument, we propose the following 
functional form for extrapolating the heavy meson decay 
constant from unphysical region to the physical limit: 



where 



/^(/l?^+/l>n/lf^) + 6i.. 



(35) 



(36) 



and b'p is a parameter. 

It is convenient to define the following new parameters: 




ap ~ 

bp — ^mpfp^ + b'p. 



(37) 



So we have the following explicit functional form for ex- 
trapolating the lattice data: 



where 



= apa^ + bp, 



d^k 1 



(38) 



f2 / (2^)3 TWfc 

f d^k |kp 



(39) 



from Eqs.(27)(32), and ap and bp are the parameters to 
be determined by fitting the lattice data. 

As mentioned in Introduction, for extrapolating the 
lattice data which are outside the power-counting regime, 
the FRR technique can be used for the resummation of the 
chiral expansion. In FRR, a finite-range cutoff in the mo- 
mentum integrals of the pion loop diagrams is introduced. 
We choose two different approaches for evaluating the in- 
tegrals in Eq.(39). One is the sharp-cutoff form, ^^(yl — |fc|), 
and the other is the dipole form, yl^/(yl'^ + |fcp)^. Both of 
them make the integrals converge. When the pion mass is 
greater than the cutoff A, which characterizes the finite 
size of the source of the pion, the Compton wavelength of 
the pion is smaller than that of the source and pion loop 
contributions are suppressed as powers of yl/m^. Obvi- 
ously the dipole form is more realistic. Since the leading 
nonanalytic contribution of pion loops is only associated 
with the infrared behavior of the integrals in Eq.(39), it 
does not depend on the details of A. 

Since the lattice simulations are performed on a finite 
volume grid, the finite size effects should be taken into 
account [15|. In the finite periodic volume, the available 
momenta k are discrete: 



27rn 
aL ' 



(40) 
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where L is the number of lattice sites in the x, y, z direc- 
tion, and the integer n is in the following range: 



Since the pion momenta on the finite lattice volume are 
discrete, we should take this into account by replacing the 
continuous integral over k in Eg. ([55)1 with a discrete sum 
over fc, 

where the discrete momenta kx^ ky, kz are given in Eq. (40) . 
From Ref. [16|. the volume of lattice simulation is 20^ x 64, 
corresponding to L = 20. The smallest momentum allowed 
on the lattice is l-njaL. 

4 Extrapolation of lattice data for 
pseudoscalar heavy meson decay constant 

Lattice gauge theory is the only quantitative tool cur- 
rently available to calculate nonperturbative phenomena 
in QCD from first principles. Quark vacuum polarization 
is the most expensive ingredient in lattice QCD simula- 
tions, especially when quark masses are very small, as in 
the case of u and d quarks. The MILC Collaboration has 
established unquenched gauge configuration, which con- 
tains three flavors of light sea quarks to include the ef- 
fects of realistic quark vacuum polarization [2]. This stag- 
gered quark discretization of QCD has several advantages, 
which are offset by the fact that staggered quarks always 
come in groups of four identical flavors. In Rcfs.[3][5]|17j. 
the Kogut-Susskind action with improved flavor and ro- 
tational symmetry suitable for dynamical fermion simula- 
tions is constructed. At tree level, the action has no cou- 
plings of quarks to gluons with a transverse momentum 
component 7r/a. So the flavor symmetry violating terms in 
the action are completely removed at tree level. Then the 
rotational symmetry is improved by introducing the Naik 
term. Finally, an extra five link staple is introduced in or- 
der to cancel errors of Oia^g^) from the fattening. The 
result action is called Asq action. It can be further im- 
proved by tadpole improvement, leading to Asqtad, which 
is an order 0{a*, a?g^) accurate fermion action. 

The authors of Ref. [I] studied the i3-meson decay con- 
stant employing the MILC Collaboration unquenched gauge 
configuration, and light sea quarks are represented by the 
highly improved staggered quark action. The good chiral 
properties of the latter action allow for a much smoother 
chiral extrapolation to the physical region. The b quark 
inside B meson is treated only as the valence quark in 
simulations because its effect in the sea should be sup- 
pressed by inverse powers of the h quark mass [18] . NRQCD 
lattice action was used to deal with the valence b quark, 
which has been developed over many years [1] [12] [10] . The 
b quark is non-relativistic inside its bound state, hence a 
non-relativistic expansion of the QCD action is appropri- 
ate which accurately handles scales of the order of typical 
momenta and kinetic energies inside these states. 



Some lattice data of are obtained in Ref. 1^ from 
full or partially quenched QCD. The coarse lattice spac- 
ing a is around 0.1 2fm. In the case of full QCD, dynamical 
light quark masses are chosen as mq/ms = 0.125, 0.175, 
0.25 and 0.5, where is the strange quark mass, corre- 
sponding to <Pb = 0.516(5)(15), 0.519(5)(15), 0.517(8)(15) 
and 0.540(5)(15) (in unit of GeV^/^), respectively (the er- 
rors in the first parentheses are statistical and the second 
come from lattice spacing uncertainties). In the case of 
partially quenched QCD, dynamical light quark masses 
are chosen as mq/ms = 0.125 and 0.5, corresponding to 

= 0.506(5)(14), and 0.547(5)(15) (in unit of GeV^/^), 
respectively. Two values of in the case of fine lattices 
with a being around 0.087fm have also been accumulated, 
where the staggered valence light propagators created by 
the Fermilab Collaboration were used to treat the light 
quarks. 

We choose to extrapolate the coarse data (full and par- 
tially quenched) since the number of the fine data is only 
two. If we used the fine data, the uncertainty of the fitted 
results due to the error in the lattice data would be too 
large. Explicit lattice simulations show that over the rang 
of mass of interest to us, where the pion mass is not con- 
strained by the chiral limit, is proportional to TO.j|21j. 

In our model, there are two parameters, ap and bp, 
to be fixed in Eq. (|38)) . These parameters are related to g, 
A2 and A. g and A2 represent the interaction at the quark 
level and cannot be determined from the chiral pertur- 
bation theory for heavy mesons. From the experimental 
data for the decay width for D D*tt, we have the vari- 
ation of g between 0.51 and 0.67[22]. From Eq.(16), A2 
is related to the mass splitting between a heavy vector 
meson and a heavy pseudoscalar meson. Based on the ex- 
perimental data for B mesons, the value of A2 should be 
around —O.OSGeV^ . To see the dependence on A2 in our 
fits, we let A2 vary between -OMGeV'^ and -OmGeV'^. 
We let A vary between OAGeV and 0.6GeV, considering 
the scale of the finite size of the pion source [10] [11]. The 
two parameters, ap and bp, are determined by the least 
squares fitting method in our fit. 

We first extrapolate the full QCD data. In this case, 
the uncertainties of the parameters as and bs arc around 
90% and 10%, respectively, due to the errors in the lattice 
data. The large uncertainty of as, associated with the 
chiral correction term, is due to the lack of the lattice data 
near the physical region. In Table 1, the fitted results in 
Columns 3 and 5 (Columns 2 and 4) correspond to the 
case where the finite lattice volume effects are (not) taken 
into account. We can see that the difference between these 
two cases is not very large. 

The fitted result for fs is obtained from Column 5 
in Table 1, where the dipole form factor is used and the 
correction from the discretization of the pion momentum 
is considered: 

/i3 = 214.2 ± 13.3 ±0.8Afey, (43) 

where the first uncertainty is from the errors in the lattice 
data, and the second is from the variation of parameters 
(5, A2, A) in our model. 
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If we use the sharp cutoff scheme, the fitted result for 
fs is 

/s = 215.1 ±15.8 ±4.1MeF. (44) 

We then extrapolate the full QCD data together with 
the two partially quenched data (although, in principle, 
one should use quenched chiral Lagrangian to extrapolate 
these two data). As expected, we find that the uncertain- 
ties which are due to the errors in the lattice data for ag , 
bg, and the extrapolated _B-meson decay constant are all 
reduced since the number of the lattice data increases (see 
Table 2). 

The fitted result for /b is obtained from Column 5 
in Table 2, where the dipolc form factor is used and the 
correction from the discretization of the pion momentum 
is considered as before: 

fs = 209.4 ± 9.7 ± l.OMeV, (45) 

where the first uncertainty is from the errors in the lattice 
data, and the second is from the variation of parameters 
(g, A2, A) in our model. 

If we use the sharp cutoff scheme, the fitted result for 
fs is 

/s = 210.5 ± 11.4 ±5. 7Mey. (46) 

It can be seen that the extrapolated results for the B- 
meson decay constant are consistent with each other in 
the above two extrapolation cases and in the two cutoff 
schemes if we consider the uncertainties of the extrapo- 
lated results. We can also see that the uncertainties due to 
variations of our model parameters in the dipole scheme, 
which is more realistic, are smaller than those in the sharp 
cutoff scheme. 



5 Summary and discussion 

It is known that the strong interactions are constrained 
by the heavy quark symmetry when the heavy quark mass 
goes to infinity, and by the chiral symmetry when the light 
quark mass approaches zero. The chiral perturbation the- 
ory can be applied when the light quark mass (or the pion 
mass) is small enough. In the region where the light quark 
mass is bigger than about 60 MeV, the non-relativistic 
constituent quark model is appropriate. Base on these, we 
propose a phenomenological functional form to extrapo- 
late the lattice data for the B-meson decay constant to the 
physical region, combining the chiral perturbation theory 
and the non-relativistic constituent quark model. 

We evaluate pion loop contributions when is small 
with the aid the chiral perturbation theory for heavy mesons. 
This leads to correct nonanalytic chiral behavior of fs in 
the chiral limit. Since the lattice data for the B-meson de- 
cay constant is outside the power-counting regime, instead 
of simply using the chiral logarithms as done in Ref.fl], 
we employ the finite-range regularization technique in or- 
der to resum higher order terms of the chiral expansion. 
At large pion masses, we introduce a constant parameter 
to fit the lattice data for <Pb (which appear in the re- 
gion where the constituent quark model is suitable to be 



applied) based on the non-relativistic constituent quark 
model. This is also different from the term linear in light 
quark mass used in Ref.[lJ. It has been pointed out that 
in large pion mass region, the extrapolation based on the 
constituent quark model is more reasonable than the sim- 
ple linear extrapolation[TT]. The finite lattice volume ef- 
fects are taken into account by replacing the continuous 
integral with a discrete sum. In our fit, the two parameters 
are determined by the least squares fitting method. 

In Ref . [1] , the lattice data for the i3-meson decay con- 
stant (both full QCD and partially quenched data) are 
extrapolated based on the staggered chiral perturbation 
theory with power counting being used. The result of the 
extrapolated B decay constant at the physical pion mass 
is 

/b = 0.216(9)(19)(4)(6)Gey, (47) 

where the errors, from left to right, are statistics plus 
scale plus chiral extrapolations, higher order matching, 
discretization, and relativistic corrections plus tuning 
respectively. 

It can be seen that the central values for the extrapo- 
lated i3-meson decay constant in Eqs.([45 |) (|46 p are smaller 
than that in Eq. (|T7)) . However, because of the large uncer- 
tainties due to the errors in the lattice data and the model 
parameters etc. the extrapolated _B-meson decay constant 
in our work is still consistent with that in Ref.yy. 

Now we discuss the uncertainties in our model. We 
have three parameters, A2, g and A, where the first two 
are related to the color-magnetic-moment operator at or- 
der 1/toq in HQET and the interaction between heavy 
mesons and Goldstones in chiral perturbation theory, re- 
spectively. Appropriate ranges for them (obtained from 
the comparison with the experimental data) are used in 
our fit process. We did not take into account the O(a^) 
contributions in the Kogut-Susskind action, which could 
give 7% corrections to our result [I]. In our approach, sea 
quark mass is not extrapolated. Usually quenched lattice 
QCD gives about 90% contribution to physical quantities, 
hence sea quark contribution is at about 10% level. There- 
fore, one expects the uncertainty from sea quark extrap- 
olation is about 1 ~ 2%. Taking these into account, our 
extrapolated result in Eqs. (j43ll46|) may have another 8% 
uncertainty. In our future work we will investigate these 
issues carefully in order to reduce these uncertainties. 
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Table 1. Fitting parameters, the fitted results in the sharp-cutoff and the dipole schemes and the (relative) uncertainties caused 
by the errors in the lattice data in the full QCD case are shown. The variations of the values in the second to the fifth columns 
are due to the variations of the parameters in our model. 



Form factor 


sharp cutoti 




dipole 




Volume 


infinite 


finite 


infinite 


finite 


as (GcV"/") 


0.3644-0.9395 


0.3759-1.0388 


0.5068-1.1943 


0.5098-1.1965 


Aflis/as (%) 


89.30-91.38 


89.36-91.62 


88.72-89.83 


88.75-89.84 


bB (GeV3/2) 


0.5715-0.5928 


0.5712-0.5910 


0.6017-0.6368 


0.5981-0.6275 


Abg/bB (%) 


7.86-10.58 


7.69-1.05 


11.80-15.88 


11.34-14.81 


$ (GeV^/^) 


0.4944-0.4994 


0.4844-0.5035 


0.4971-0.5000 


0.4900-0.4936 


(GeV^/^) 


0.0227-0.0275 


0.0195-0.0364 


0.0221-0.0248 


0.0277-0.0305 


A^/$ (%) 


4.54-5.57 


3.87-7.52 


4.42-4.99 


5.61-6.23 


fs (GeV) 


0.2153-0.2175 


0.2110-0.2193 


0.2165-0.2178 


0.2134-0.215 


A/b (GeV) 


0.0099-0.0120 


0.0085-0.0159 


0.0096-0.0108 


0.0120-0.0133 


A/b//b (%) 


4.54-5.57 


3.87-7.52 


4.42-4.99 


5.61-6.23 



Table 2. Fitting parameters, the fitted results in the sharp-cutoff and the dipole schemes and the (relative) uncertainties caused 
by the errors in the lattice data in the case where both the full QCD and partially quenched QCD data are fitted are shown. 
The variations of the values in the second t(i the fifth cdluniiis are due to variation of tin; jjaraiiit^tcn's in our model. 



Form factor 


sharp cutoff 


sharp cutoff 


dipole 


dipole 


Volume 


infinite 


finite 


infinite 


finite 


as (GeV«/^) 


0.5160-1.3369 


0.5075-1.4148 


0.6810-1.6116 


0.6851-1.6146 


Aos/os (%) 


49.47-52.91 


47.35-47.78 


47.23-47.43 


47.24-47.44 


bB (GeV3/2) 


0.5944-0.6238 


0.5878-0.6160 


0.6295-0.6758 


0.6247-0.6633 


A6b/6b (%) 


6.22-8.69 


5.31-7.16 


8.03-10.66 


7.72-9.98 


^B (GeV^/^) 


0.4839-0.4910 


0.4707-0.4968 


0.4971-0.5000 


0.4788-0.4835 


A#s (GeV3/2) 


0.0177-0.0217 


0.0144-0.0262 


0.0163-0.0182 


0.0202-0.0223 


A<Pb/<Pb (%) 


3.62-4.48 


2.90-5.58 


3.32-3.72 


4.18-4.66 


fB (GeV) 


0.2105-0.2137 


0.2048-0.2162 


0.2163-0.2176 


0.2084-0.2104 


A/s (GeV) 


0.0077-0.0094 


0.0063-0.0114 


0.0071-0.0079 


0.0088-0.0097 


AfB/fB (%) 


3.62-4.48 


2.90-5.58 


3.32-3.72 


4.18-4.66 
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Fig. 1. Heavy meson decay constant to one pion loop order. The black squares represent the weak current in Eq.(22), and 
dashed lines denote the pion. 



P P 

(a) 



(c) 



P P 

(d) 



10 



X.-H. Guo, M.-H. Weng: Chiral extrapolation of lattice data for _B-meson decay constant 



Fig. 2. Extrapolation of the full QCD lattice data for the B-meson decay constant. Dashed line (dotted line) corresponds to 
the sharp cutoff scheme and infinite (finite) volume; dash-dot line (solid line) corresponds to the dipole scheme and infinite 
(finite) volume. The vertical solid line corresponds to the physical pion mass. 
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Fig. 3. Extrapolation of both the full QCD and partially quenched QCD lattice data (which are represented by squares and 

diamonds respectively) for the B-mcson decay constant. Dashed line (dotted line) corresponds to the sharp cutoff scheme and 
infinite (finite) volume; dash-dot line (solid line) corresponds to the dipole scheme and infinite (finite) volume. The vertical solid 
line corresponds to the physical pion mass. 
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